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ABSTRACT 

M 87 is a nearby radio galaxy that is detected at energies ranging from radio to 
very-high-energy (VHE) gamma-rays. Its proximity and its jet, misaligned from our 
line-of-sight, enable detailed morphological studies and extensive modeling at radio, 
optical, and X-ray energies. Flaring activity was observed at all energies, and multi- 
wavelength correlations would help clarify the origin of the VHE emission. In this paper, 
we describe a detailed temporal and spectral analysis of the VERITAS VHE gamma- 
ray observations of M87 in 2008 and 2009. In the 2008 observing season, VERITAS 
detected an excess with a statistical significance of 7.2 standard deviations (<r) from M 87 
during a joint multi-wavelength monitoring campaign conducted by three major VHE 
experiments along with the Chandra X-ray Observatory. In February 2008, VERITAS 
observed a VHE flare from M87 occurring over a 4-day timespan. The peak nightly 
flux above 250 GeV was (1.14 ± 0.26) x 10 ^crtT^s , which corresponded to 7.7% of 
the Crab Nebula flux. M 87 was marginally detected before this 4-day flare period, and 
was not detected afterwards. Spectral analysis of the VERITAS observations showed 
no significant change in the photon index between the flare and pre-flare states. Shortly 
after the VHE flare seen by VERITAS, the Chandra X-ray Observatory detected the 
flux from the core of M 87 at a historical maximu m, while the flux from the nearby knot 



HST-1 remained quiescent. lAcciari et al.l (|2009i ) presented the 2008 contemporaneous 



VHE gamma-ray, Chandra X-ray, and VLBA radio observations which suggest the core 
as the most likely source of VHE emission, in contrast to the 2005 VHE flare that was 
simultaneous with an X-ray flare in the HST-1 knot. In 2009, VERITAS continued its 
monitoring of M87 and marginally detected a 4.2 a excess corresponding to a flux of 
~ 1% of the Crab Nebula. No VHE flaring activity was observed in 2009. 

Subject headings: gamma rays: galaxies — galaxies: individual (M 87, VER J1230+123) 



1. Introduction 

M87 is an FR I radio galaxy located at a distance of 16.7 Mpc, near the center of the Virgo 
cluster. It has been observed at all wavelengths ranging from radio to VHE gamma rays. Its core is 
an active galactic nucleus (A GN) powered by a super ma ssive black hole of m ass (6.0±0.5)xl0 9 M 



(jGebhardt &: Thomaal2009i ) (see supporting material of lAcciari et al.l (120091 ) for mass corrected for 



our distance assumption), which is the source from which the first plasma jet emission was observed 



N Now at Institut fur Physik und Astronomie, Universitat Potsdam, 14476 Potsdam-Golm, Germany; DESY, Pla- 
tanenallee 6, 15738 Zeuthen, Germany 

§ Now at Los Alamos National Laboratory, MS H803, Los Alamos, NM 87545 
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(|Curtislll918r ). Most of the known extragalactic VHE sources are blazars (|TeVCatJl2007l ). AGN with 
a jet aligned close to the line-of-sight; on the cont rary, the jet of M8 7 is misalign ed. Apparent 



super luminal motion was observed in both the radio (jCheung et al.ll2007l ) and optical (jBiretta et al 



1999 ) bands for different features along the jet, constraining the jet orientation to less than 30° 
from the line-of-s ight at the location of the HST-1 knot, which is located 0.86" from the core 
(|Harris et alJboOfll ). 



VHE emission from M 87 was first reported by t he HEGRA collaborati on at a statistical signifi- 
cance of 4.1 a during their 1998 — 1999 observation s (jAharonian et al.ll2003l ). This was confirmed by 
the H.E.S.S. collaboration (jAharonian et al.ll2006l ). which additionally reported year-scale flux vari- 
ability. The observed variability provides a size constraint on the VHE emission production region 
and disfavors larg e-scale gamma-ray production models, such as t he dark matter annihilation model 
( Balt z et al hood ) and the interacting cosmic-ray proton scenario ( Ipfrommer fc Ensslinlhoosl ) which 
predict steady gamma-ray emission. However, the angular resolution of ima ging atmospheric 



Chere nkov telescopes (IACTs) is insufficient to resolve any structure in M87. lAharonian et al 
(|2006l ) also reported fast (2-day scale) variability during a high state of gamma-ray activity in 
2005, which further constrains the VHE emission size and favors the immediate vicinity of the M87 
black hole as the VHE production site; on the other hand, Chandra X-ray observations at the time 
of the VHE flare indicated a different scenario. Strong flux variability from both the core and 
iST-1 in the energ y range 0.2 — 6keV has been detected by Chandra M87 monitori ng since 2002 



Har ris et al.ll2003l ). At the same period of the flare observed by H.E.S.S. in 2005, lHarris et al 



2003) reported an X-ray flux from HST-1 at more than 50 times the intensity observed in 2000, 



along with flaring activity observed in ultraviolet and radio wavelengths, suggesting the HST-1 
knot as a more likely source of VHE emission than the core. 

The All-Sky Monitor (ASM) on the Rossi X-ray Timing Explorer (RXTE) has provided daily 
monitoring of M87 in the energy range 2 — 12keV since 1996. However, ASM cannot resolve the 
core and the HST-1 knot. A year-scale correlation between the annual A SM /RXTE X-ray flu x and 
the VHE gamma-ray flux recorded over the past 10 years was reported by lAcciari et al.l (|2008l ). The 
core was presented as the more favorable VHE production site over HST-1 based on this year-scale 
correlation, and the non-detection of the 2005 flare in the HST-1 knot by ASM. However, no year- 
scale correlation has been seen with Chandra data of the M 87 core or the HST-1 knot over the past 
5 years in the 0.2 — 6keV range. The ASM/RXTE quick- look results do not show any correlated 
activity with VHE observations at shorter time scales, as a result of the limited sensitivity of ASM. 

The proximity of M 87 a nd its misaligned jet e n able high resolution stud ies of its jet structures 
in radio, optical, and X-rays ([Marshall et al.ll2002l ) (jPerlman fc Wilson! 120051 ). The jet morphology 
is similar in those wavebands, an d leptonic synchrotron radiation is favored as the process for non- 
thermal emission within the jet (IWilson fc Yangl 120021 ) . Based on the multiwavelength correlated 
activities and variabili ty studies, the favore d candidate for VHE emission is the small-scale (sub- 
arcsecond) jet region. iReimer et al.1 ([2004]) propose a synchrotron-proton blazar model in which 
protons are accelerated to energies above EeV and emit gamma rays via muon/pion synchrotron 
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or proton synchrotron radiation. However, this model requires a strong magnetic field to acceler- 
ate the protons to such high energies, and the predicted VHE spectrum is steeper than observed. 
Several leptonic models involving synchrotron and inverse Compton (IC) radiation have also been 
suggested, with the multi-co mponent emission originating in the inner jet: there is the model by 
Georganopoulos et al.l ((2005) in which energetic electrons IC scatter off o f sync hrotron photons 
produced downstream in the decelerating jet; a scenario by iLenain et al.l (|2008l ) in which VHE 
emission is produced via the synchrotron self-Compton (SSC) process inside several similar ho- 
m ogenous compact components w hich contain more energetic electrons than the jet; and the model 
by iTavecchio fc Ghisellinil (|2008l ) , in which a fast-moving spine is surrounded by a slower- moving 
sheath, producing VHE emission via exter nal-Compton scattering. Single-component SSC emission 
is recently modelled bv lAbdo et al.l (|2009l ) with VLBA radio, Chandra X-ray, and Fermi LAT 2009 
data. However, comparing the model to archival non-flaring VHE data, it appears to underpredict 
the VHE emission. 



The vicinity of the black hole (the core) has been suggested by lLevinsonl (|2000l ) and lNeronov &; Aharonian 



(|2007l ) in the black hole magnetosphere model, in which gamma-ray photons are produced by elec- 
trons accelerated by the electromagnetic field of the black hole. The HST-1 knot, 0.86" away from 
the core, has also been demonstrated as a possibl e location for jet rec onfinement where photons 
can be upscattered to TeV energies via IC process ( Stawarz et al.ll2006l ). 



In 2008, the VHE gamma-ray experiments H.E.S.S., MAGIC, and V ERITAS took part i n 
a joint multi- wavelength monitoring campaign of M87 along with Chandra (jBeilicke et al.ll2008l ). 
The VHE observations were closely coordinated in order to guarantee a reasonable coverage around 
the Chandra pointings and a well sampled VHE light curve during the first half of 2008. During 
this joint monitoring campaign, MAGIC reported flaring activities during a 13-day observation 
period between January 30 (MJD 54495) and Fe bruary 11 (MJD 5 4507), with day-scale variability 
occurring throughout the duration of the flare ([Albert et al.l 120081 ) . Subsequently, the VERITAS 
collaboration triggered intensified observations of M 87, and detected a 4-day flare from February 9 
(MJD 54505) to February 13 (MJD 54509). The VHE and Chandra X-ray li ght curves of the join t 
campaign, along with a coincident VLBA radio light curve, are presented in lAcciari et al.l ( 20091 ) . 
which also include a discussion on the radio/ VHE gamma-ray correlation as evidence that the VHE 
emission originates from the core. 

The VHE flux monitoring campaign continued in 2009 with MAGIC and VERITAS. In this 
paper, we present the results from two seasons (2008 — 2009) of VERITAS observations, along with 
a full analysis on the time scale of the flux variability, and a search for spectral variability in the 
2008 data set when the VHE flaring activity was observed. 
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2. Observations and Analysis 



VERITAS, the Very Energetic Radiation Imaging Telescope Array System, is an array of 
four 12 m diameter IACTs located at the Fred Lawrence Whipple Observatory on Mount Hopkins 
(31°40'N, 110°57 / W) at an altitude of 1.3 km above sea level. Each telescope has a total area of 
110 m 2 and an f/D ratio of 1.0. Each telescope camera is equipped with 499 photomultiplier tubes, 
arranged in a hexagonal lattice covering a field of view (FOV) of 3.5°. The array is sensitive from 
~100 GeV to more than 30 TeV, with an effective area of up to 10 5 m 2 and an angular resolution of 
0.1° at 1 TeV (68% containment). VERITAS can detect a source with 1% Crab Nebula flux in less 
than 50 hours, and a source with 5% Crab Nebula flux in ~2.5 hours. For more technical details 



of VERITAS, see lHolder et al.1 (|2006h 



M87 was observed with VERITAS for over 43 hours between December 2007 and May 2008 
and over 25 hours between January and May 2009, at a range of zenith angles from 19° to 41°. All 
of 2008 data and 81% of 2009 data used were taken with 4 telescopes, while 19% of 2009 data were 
taken with only a 3-telescope array. After eliminating observations in poor weather and those with 
an unstable trigger rate, 37 hours and 19 hours of good quality live time remain in 2008 and 2009, 
respectively. 

Shower images from all working telescopes are first corrected in gain a nd timing using param- 
eters obtained from the nightly laser calibration data (jHanna et al.1 120071 ) . The images are then 
passed through a two-step cleaning process that retains pixels with a signal that is several time s 



higher than the night-sky background level. Each shower image is then para metrized ( Hillas 



and the shower direction is reconstructed using the stereoscopic technique (jHofmann et al 



1985) 



19991 ) 



Events are then selected as gamma-ray- like if a t least three came ra images pass the standard cuts 
optimized for a 10% Crab Nebula flux source (jColin et al.ll2008l ). All the observations were per- 
formed in "wobble mode" where M 87 was tracked with a 0.5° offset relative to the camera such 
that the camera's FOV contain both the source region and regions for background estimation. The 
on-source region is defined by a 0.15° radius circle centered on the M87 core. All gamma-ray-like 
events within this region are summed (ON count) and the background estimated from seven identi- 
cally si zed regions reflect ed from the source region around the camera center is als o summed (OFF 
count) (jBerge et ajj[2007| ). The ON and OFF counts are then used in Formula 17 of lLi fc Mai (119831 ) 
to calculate the significance of the excess. As a standard procedure in VERI TAS, the results were 
co nfirmed by at least one independent analysis package ( 
in lAcciari et al.1 (hood ). 



Daniel et al.l 120071 ) which was presented 



-7- 



3. Results 
3.1. 2008 

The VERITAS observations in 2008 resulted in 450 excess events, corresponding to a statistical 
significance of 7.2 a. The average flux above 250 GeV is (2.74±0.61) x 10~ 12 cm -2 s" 1 , corresponding 
to 1.8% of the Crab Nebula flux. The differential energy spectrum of M87 measured by VERITAS 
in 2008 is consistent with a power- law distribution dN/dE = $ (£/TeV)- r , with $ = (5.17 ± 
0.91 stat ± im syst ) x 10" 13 cm" 2 s" 1 TeV" 1 and Y = 2.49 ± 0.19 stat ± 0.20 s2/st . The x 2 /dof for the 
powe r-law fit is 3.1/4 . 0. Th e measured photon i ndex i s consistent with re ported measurements 



from lAharonian et al.l (|2006l ). lAcciari et al.l (|2008l ). and lAlbert et al.l (|2008l ) ranging from 2.22 to 



2.62 with statistical uncertainties between 0.11 and 0.35. Table [T] lists the differential flux points 
measured by VERITAS in 2008. 

The 2008 VERITAS light curve is shown in Fig. [I] assuming the average fitted photon index 
of 2.50. The x 2 /dof of a constant rate fit to the entire dataset is 52/26, rejecting the constant 
flux hypothesis at 99.9% confidence level. On Feb 9 (MJD 54505), M87 was detected at over 4 a 
after 2 hours of observation and a trigger alert was sent out to other collaborations for intensified 
observations. The peak flux of the VERITAS dataset occured on the night of Feb 12 (MJD 54508) 
at (1.14 ± 0.26) x lO^cm^s- 1 above 250 GeV, equivalent to 7.7% of the Crab Nebula flux. The 
flare period is defined by the nights removed from the whole dataset such that the constant rate 
fit to the rest of the lightcurve reaches a x 2 /dof close to 1. The nights that meet the criteria for 
the flare period are Feb 9, 10, 12, and 13 (MJD 54505 - 54509). The dataset is then further split 
into pre-flare and post-flare period. 

During the flare period, M87 was observed for 5.4 hours of live time and was detected at 7.4 a 
(pre-trials). The flux above 250 GeV during the flare period was (7.59 ± 1.11) x 10~ 12 cm" 2 s , 
corresponding to 5.1% of the Crab Nebula flux. The 99% confidence interval for the flux during 

Table 1. M87 differential flux measured between December 2007 and May 2008 by VERITAS. 



Energy 


Differential Flux 




Significance 


[TeV] 


[cm" 2 s- 1 TeV" 1 ] 




w 


0.24 


(1.99 ±0.54) x 10" 


ii 


3.68 


0.42 


(3.56 ± 1.16) x 10" 


12 


3.07 


0.75 


(1.51 ±0.34) x 10- 


12 


4.44 


1.33 


(1.70 ±0.95) x 10- 


13 


1.79 


2.37 


(6.02 ± 3.40) x 10- 


14 


1.78 


4.22 


(1.72 ± 1.13) x lQ- 


11 


1.52 
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Fig. 1. — VERITAS night-by-night VHE gamma-ray flux and Chandra X-ray flux ([Harris et al 



2009|) from the core and from the HST-1 knot of M 87 in 2008. The flare period (MJD 54505-54509, 
Feb 9 — 13) seen in VHE gamma rays coincides with a historically high state of the core in X- 
rays, while the HST-1 knot remains in a low flux state. The black slant-lined area shows the 99% 
confidence intervals of the pre-flare period; the black horizontal-lined area shows the 99% confidence 
intervals of the flare period; the black line wi th the arrow indicates the upper limit of the post-flare 
period at 99% confidence level (|Helendll983l ). 
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the flare period (see horizontal-lined area between MJD 54505 and 54509 in Fig. [I]) is between 
4.73 x 10~ 12 cm -2 s -1 and 10.45 x 10~ 12 cm" 2 s" 1 



Before the flare period (MJ D 54448 — 54503), w hich included two nights immediately after the 



major flare observed by MAGIC (lAlbert et al.ll2008l ). M87 was marginally detected with 13.8 hours 



of live time at a flux of (3.30±0.68) x 10~ 12 cm" 2 s" 1 above 250 GeV. The 99% confidence interval for 
the flux during the pre-flare period (see slant-lined area between MJD 54448 and 54503 in Fig. [I]) is 
between 1.55 x 10 -12 cm~ 2 s _1 and 5.05 x 10 -12 cm -2 s _1 . M87 was observed for an additional 17.2 
hours of live time after t he flare perio d and was not detected. An upper limit of 2.1 x 10 -12 cm -2 s -1 
at 99% confidence level (jHelendl 19831 ) is established for the post-flare period, corresponding to 1.4% 
of the Crab Nebula flux. Assuming a normal distribution for the flux measurement, the flare flux 
is higher than the pre-flare flux at 99.95% confidence level, and the pre-flare flux is higher than the 
post-flare flux at 99.8% confidence level. All of the above calculations were performed under the 
assumption of a constant photon index of 2.50. 

The spectra from different periods are plotted in Fig. [2] and there is no indication of a spectral 
cut-off. The power-law fit parameters from the overall, flare, and pre-flare spectra are displayed in 
Table [2) No significant difference in photon index and photon energy distribution is found between 
the flare and pre-flare states. 

To test if there is any correlation between the spectral flux and photon index measured during 
different activity levels corresponding to the c ore, a linear test functio n of the form r = pn + Pi x 



loglO $o was fitted to past measurements by lAharonian et al.l (120061) in 2004 and lAcciari et al 



(|2008l ) in 2007, and the high/low states measured bv klbert et al.l ifcood ) and by VERITAS in 2008 
(see Fig. [3]). The measurement by HESS in 2005 is excluded in the fit due to possible contamination 
from the HST-1 flare. The x 2 /dof of the linear fit is 1.7/4 and the corresponding probability of a 
correlation between photon index and flux is 78.6%. The fit parameter p\ is consistent with zero. 
Fitting a constant photon index gave a \ 2 /dof of 2.7/5 and a corresponding probability of 74.6%. 
No significant correlation is found between the photon index and the flux. 

At the time of the flare detected by VERITAS, Chandra measured the core X-ray flux at a 
historical maximum, at 4.1 a above the mean core flux between 2000 and 2009, and exceeding the 



Table 2. Power-law fit parameters of the VHE spectrum of M87 in different epochs during 2008. 



Data Set 




r 


X 2 /dof 




[10- 13 cm- 2 s- 1 TeV- 1 ] 






2008 overall 


5.2 ± 0.9 


2.5 ± 0.2 


3.1/4 


flare state 


15.9 ± 2.9 


2.4 ± 0.2 


3.9/4 


pre-flare state 


5.6 ± 1.5 


2.5 ± 0.3 


1.3/4 
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Fig. 2. — Spectra of M87 measured at different periods by VERITAS. Data points less than 1.5 a 
are represented as upper limits. The power-law fit parameters are listed in Table [21 The photon 
energy distribution is consistent throughout different activity levels, and there is no indication of a 
spectral cut-off. 
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Fig. 3. — Photon in dex versus differential flux as reported by HESS in 2004 and 200 5 
(|Aharonian et al.1 bood ). b y MAGIC in 2008 (overall, flare, and non-flare) (|Albert et allboosl ). 
and by VERITAS in 2007 (jAcciari et alJbood ) and 2008 (overall, flare, and pre-flare). All spectra 
are compatible within their statistical errors. A linear coorelation fit between photon index and 
flux has the same probability as a constant photon index fit. 



flux from the HST-1 knot (see Fig.[T]and|3D ([Harris et al.l l2009). The peak core X-ray flux measured 
by Chandra during the observation period of VERITAS was 2.3 times the average core X-ray flux, 
with the average calculated excluding this flare data point. While the HST-1 knot X-ray flux during 
this period fluctuated ±10% from the average and is relatively steady when c ompared to the cor e 
emission. Fig. [5] shows the fractional change per year (fpy) first presented in lHarris et al.l ([2009). 
Around the 2008 VHE gamma-ray flare period, the Chandra core fpy is more v ariable than that o f 
the Chandra HST-1 fpy measurement. The definition of fpy is repeated below ([Harris et al.ll2009l ). 



fpy 
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Fig. 4. — M 87 y early VHE gamm a-ray (jAharonian et al.l l2003l) (lAharonian et alj [2006) 



(jAcciari et all l2008h (| Albert et al.ll2008l ) and X-ray fluxes (Harris et al.ll2009h . The VHE gamma- 
ray flux is for energy >730GeV due to the original flu x scale used in the HEGRA paper. This 
yearly flux plot is first presented in (jAcciari et al.ll2008l ) and is now updated with 2008 and 2009 
data. Grey areas represent the range of variable VHE fluxes observed that year to give a more 
accurate picture of the flux level of M87. 
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where / is the X-ray intensity measured by Chandra and At is the time between the two intensity 
measurements in units of year. If the X-ray intensity is increasing (i.e. I2 > Ii), then i = 1, j = 2; 
if the X-ray intensity is decreasing (I2 < then i = 2,j = 1. 



3.2. 2009 

The VERITAS observations in 2009 resulted in 134 excess events at a significance of 4.2 a. The 
X 2 /dof of a constant rate fit to the entire dataset is 23/18 and no significant variability was observed 
in the 2009 data set. The average flux above 250 GeV is (1.59 ± 0.39) x 10~ 12 cm -2 s _1 assuming a 
photon index of 2.50. This corr esponds to 1-1% o f the Crab Nebula flux and is consistent with the 



reported Fermi-LAT spectrum (jAbdo et al.ll2009l ). The 2009 flux is below the 2008 pre-flare flux at 



the 98.5% confidence level, but above the 2008 post-flare flux at the 87.6% confidence level. Fig.[5] 
shows the light curves of both 2008 and 2009 observed by VERITAS and the X-ray flux measured 
by Chandra for the core and the knot HST-1, along with the fpy of Chandra X-ray flux. 

During the 2009 observation period, Chandra took a flux measurement every 40 — 50 days. The 
HST-1 flux measured by Chandra is steadily declining at a rate of 5 — 10% when compared to the 
previous flux (51 The core flux, however, is mostly increasing by as much as 18%. The lack of 
VHE activity during this period suggests a possible association between VHE flares and significant 
changes in the X-ray flux such as the one seen in 2008 where the core fpy is more variable than in 
2009 (see Fig.©. 



4. Discussion 



The proximity of M87 and its misaligned jet have enabled the study of its jet morphology 
in a broad range of energies. Flaring activities from indiv idual jet features have been observed in 
radio, optical, and X-rays in parallel (jCheung et al.1 120071 ). Modeling of the particle acceleration 
yields several possible VHE emission origins. Even though the VHE gamma-ray techn i que ca nnot 
resolve indiv i dual f eatures of M 87, the rapid variability reported by lAharonian et al.1 (|2006l ) and 
Albert et al.l (|2008l ) has constrained the size of the VHE emission region to < 2.6 5 R s where 5 
is the relativistic Doppler factor and R s the Schwarzschild radius of the M87 black hole (R s ~ 
10 15 cm). During the 2008 joint monitoring campaign of M87, VERITAS observed a gamma-ray 
flare in February 2008 which spanned 4 days, constraining the emission region size to R ^ Rvar — 
5cAt/(l + z) ~ 5 10 16 cm. Radio observations have shown evidence that c harged particles ar e 



accelerated in the immediate vicinity of the black hole closer than 100 R s (lAcciari et al.1 120091 ) . 
Combining these radio observations with the size constraint from the VERITAS 2008 data, we can 
constrain the relativistic Doppler factor S < 10 with the flaring activity measured by VERITAS. 



The VHE emission size constraint narrows down the most probable VHE emission location 
to either the unresolved core region or the HST-1 knot. X-ray time scale analysis presented by 
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Fig. 5. — Upper panel: VERITAS night-by-night VHE gamma-ray flux and Chandra X-ray flux 
from the core and from the HST-1 knot of M87 in 2008 ([Harris et al.ll2009) and 2009. Low er panel- 
Chandra X-ray flux fractional change per year, see text for definition ([Harris et al.l 12009 ). 



Harris et al.l (|2009i ) suggests the X-ray emission size of the core to be smaller than that of HST-1. 
However the upper limits on the X-ray emission region size from the time scale analysis are still 
larger than the VHE emission size inferred from the VHE variability time scales reported previously. 
Therefore, both HST-1 and the core region remain as the possible origins of VHE emission. 

In 2009 VERITAS detected M87 similar to the 2008 pre-flare level. The X-ray time scale 
analysis is repeated for the 2009 Chandra data and the fractional change per year (fpy) for the 
core in 2009 is smaller than in 2008 during the VHE gamma-ray flare. The fpy of the knot HST-1 
appeared to be negative throughout 2008 and 2009, and at a smaller amplitude than the core. This 
may be an indication of correlation between large X-ray flux changes and flaring activities in VHE 
gamma rays. It should be noted however, that large fpy in the core was observed by Chandra as 
well in 2007 with no corresponding VHE gamma-ray flare. A similar analysis is not performed for 
the 2005 HESS flare due to contamination from the knot HST-1 in Chandra data. 



The compactness of the particle accelerators operating in the vicinity of the supermassive 
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black hole and the absence of a signi ficant cut- o ff in t he spectrum imply that the p article accel- 
eration mechanism is highly efficient. iLevinsonl (|200fj ). iNeronov &: Aharonianl (|2007l ). and others 
have argued that acceleration mechanisms similar to those in pulsar magnetosphere models can also 
operate in the bl ack hole magnetosphere. Th e black hole magneto s phere model and the two-zone 
leptonic models (|Georganopoulos et al.l 120051 ) (|Lenain et al.l 12008 ) (jTavecchio k, Ghisellinil |2008| ) 
can reproduce the VHE spectrum and the flux variability well. While no constraints can be placed 
on the VHE range of these models, since VERITAS and other VHE gamma-ray instruments did not 
observe a cut-off in the spectrum of M 87, the Fermi gamma-ray telescope, sensitive from MeV to 
GeV energies, can potentially provide constraints on the magnetic field strengt h and a mea s ure o f 
low-energy gamma-ray photon flux from the synchrotron/curvature radiation. lAbdo et al.l ([2009) 
presented the first year of Fermi-LAT observations of M87 in 2009 during its quiescent state. Con- 
tinual monitoring in all wavelength is essential for the modelling of the spectral energy distribution 
of M87. Future instruments such as CTA/AGIS, with their improved sensitivity, can potentially 
detect shorter time scale variability and further constrain the size of the VHE emission region. 



Begelmanl (120091 ) provided a concise commentary on the M87 2008 multi-wavelength results 



and the importance of multi-wavelength observations to understanding particle acceleration mech- 
anisms near black holes. The high-resolution imaging of radio telescopes, combined with VHE 
gamma-ray observations from ground-based instruments, have shown the first association between 
a V HE gamma-ray flare and an increase in radio flux coming from a region very close to the black 
hole (jAcciari et al.ll2009l ) . This coincident multi- wavelength cover age led to seve r al mo deling discus- 
sions of the joint VHE and radio light curves (see supplement of lAcciari et al.1 (120091 ) ). Even with 
M87 day-scale VHE flux variability and additional constraints from radio and X-ray observations, 
there remains several plausible models, which explain how particles are accelerated to very high 
energies near the black hole and how the consequent radiation is able to reach us. Multi- wavelength 
monitoring work is being continued, and is essential to address these questions. 
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